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Abstract 

We elaborate on previous results concerning the positivity of the Euclidean path integral measure 
for low-energy modes in dense fermionic matter. We show that the sign problem usually associated 
with fermions is absent if one considers only low-energy degrees of freedom. We describe a method 
for simulating dense QCD on the lattice and give a proof using rigorous inequalities that the 
color-flavor locked (CFL) phase is the true vacuum of three flavor, massless QCD. We also discuss 
applications to electronic systems in condensed matter, such as generalized Hubbard models. 
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I. INTRODUCTION 



Euclidean quantum chromodynamics (QCD) with a non-zero chemical potential has a 
complex measure, which has made lattice simulation particularly difficult (For recent 
lattice investigation of the QCD phase boundary at finite density, see 0, This problem 
is often referred to as the sign problem, because by appropriately grouping terms quantities 
such as the partition function can be written as a sum over real, but potentially negative, 
terms. (That this grouping can be accomplished is in many systems a consequence of 
a discrete symmetry such as parity or time-reversal invariance.) Indeterminate signs are 
enough to preclude use of importance sampling, the main technique for speeding up Monte 
Carlo integration. It is important to note that while the sign problem often arises in systems 
of fermions, it is neither inevitable nor inescapable. For example, in QCD at zero chemical 
potential and in the Hubbard model at half filling one can organize the sum so that terms 
are real and positive. For recent work on the sign problem, see 0. 



Analytical work in color superconductivity [5] has demonstrated a rich phase structure 
at high density, and stimulated interest in QCD at non-zero baryon density. Several experi- 
ments have been proposed to probe matter at density of a few times nuclear matter density 
Even rudimentary information about the behavior of dense matter would be useful to 
the experimental program, as well as to the study of compact astrophysical objects such as 
neutron stars. In an earlier paper 0], we showed that QCD near a Fermi surface has positive, 
semi-definite measure. In the limit of low energies, the contribution of the remaining modes 
far from the Fermi surface can be systematically expanded, using a high density effective 
theory previously introduced by one of us jM 0| ■ This effective theory is sufficient to study 
phenomena like color superconductivity, although quantities like the equation of state may 
be largely determined by dynamics deep in the Fermi sea. 

The expansion about the Fermi surface is in powers of where /i is the chemical 
potential. For this expansion to be controlled, the ultraviolet cutoff of our effective theory 
must be less than /i, or equivalently the scale of the physics of interest must be small 
relative to the chemical potential. In QCD at asymptotic density, the superconducting gap 
is exponentially small, so this condition is satisfied. However, it is also quite possible that 
at intermediate densities (e.g., those inside a neutron star) the gap is somewhat smaller 
than fi, providing us with an additional small dimensionless parameter. Even if this is 
not the case, the power expansion of the effective theory is qualitatively different from 
the usual perturbation in as, and therefore worth exploring. Finally, we note that models 
(e.g., of electronic systems) which are not asymptotically free may exhibit strongly coupled 
quasiparticle excitations even at high density. The results described here still apply to such 
systems and may be of use in their simulation. 

This paper is organized as follows. We begin in section 2 with simple examples in (1+1) 
and higher dimensions which illustrate how the effective Fermi surface description can have 
positive measure even if the original model has a sign problem. In section 3 we review the 
results from our previous paper, and include some discussion of how they apply to electronic 
systems that may arise in condensed matter. In section 4 we describe how they can be applied 
to lattice simulations of dense matter. In section 5 we discuss correlator inequalities (also 
known as QCD inequalities) that result from positivity, and how they restrict the possible 
ground states of QCD. In section 6 we conclude with a summary and future prospects. 
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II. EXAMPLE: (1+1) DIMENSIONS AND BEYOND 



We begin with an example that illustrates the basic ideas in a simple setting. Consider 
the Euclidean (1+1) action of non-relativistic fermions interacting with a gauge field A 

S = j drdx i)l[{-dr + i(t) + tp) - e{-id^ + A)]^:^ (1) 

where e{p) is the energy as a function of momentum (e.g. e{p) ~ ^ + " " " )■ Iii © below 
we may consider it as a function of the operator {—idx + A). The dispersion relation in the 
presence of the chemical potential ep is: E{p) = e{p) — ep, and a low energy mode must 
have momentum close to ±Pf, where e(±pi?) = ep. The Fermi surface in (1+1) dimensions 
is reduced to the two points p = ±pp. Near these points we have 

E{p±pp)^±vpp , (2) 

where vp = dE/dp\pp is the Fermi velocity. 

The action (0) is not obviously positive. In fact, the operator in brackets [ ■ ■ ■ ] clearly 
has Hermitian as well as anti-Hermitian components, and hence complex eigenvalues. 

Let us assume that the gauge field has small amplitude and is slowly varying relative to 
the scale pp. We will extract the slowly varying component of the fermion field to construct 
a low energy effective theory involving quasiparticles and gauge fields. This effective theory 
will have positive, semi-definite determinant. 

First, we extract the quasiparticle modes (we suppress the spin index in what follows) 

^|Jix,T)=^PLe+'P^^ + ^rc-'^^^ , (3) 
where the functions 4'l,r are slowly varying. To simplify the action, we use the identity 

^±iPFX E{-idx + A) e^'P^"" ^{x) ^ ± vp{-idx + A)V'(x) , (4) 

to obtain [l3| 

5'cff = J drdx ip]^{—dr + iff) + idx — A)iIjl + i^iii—dr + i(f> — id^ + A)ipji . (5) 

We can write this in a more familiar form by introducing the Euclidean (1+1) gamma 
matrices 70,1,2 , which are Hermitian and can be taken as 7^ = cxj where a are the Pauli 
matrices. Using '4'l,r = |(1 ±72)"^ we obtain 

S^s = j drdx i''j^{d^ + iA^)^p = j drdx tppip . (6) 

Since the gamma matrices are Hermitian, and the operator (9^ + iA^) is anti-Hermitian, 
the operator ^ in (jH)) has purely imaginary eigenvalues. However, because 72 anticommutes 
with p, the eigenvalues come in conjugate pairs: given p 4> = X(f>, we have 

P (720) = -72-P = -72A0 = -A(720„) . 
Hence the determinant det^ = Yl A*A is real and positive semi-definite. 
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Thus, by considering only the low-energy modes near the Fermi points of the original 
model (0), we obtain an effective theory with desirable positivity properties. Note that it 
is necessary that the interactions (in this case, the background gauge field A) not couple 
strongly the low-energy modes to fast modes which are far from the Fermi points. This is 
a reasonable approximation in many physical situations, where it is the interactions among 
quasiparticles that are of primary interest. In what follows, we will apply this basic idea to 
more complex models such as QCD. 

It is straightforward to go beyond (1+1) dimensions. Consider an electron system, de- 
scribed by 

£ = V>t[z9t-e(p)]V + /iV'V, (7) 

where e(p) is the electron energy, a function of momentum p. It is interesting to note that the 
non-relativistic system already has a sign problem even at the zero density, = 0, though 
the free case does not suffer this, thanks to the separation of variables. In fact, it is quite 
unusual to have a system like vacuum QCD which has no sign problem. In Euclidean space 
the electron determinant is 

M = -a, - e{p) + /i. (8) 

The first term in operator (jHJ is anti-Hermitian, while the rest are Hermitian. Since there is 
no constant matrix P in the spin space that satisfies M'^ = P M P~^, it has a sign problem 
in general. 

Let us decompose the fermion momentum as 

P = Pf + L (9) 

Again, the Fermi momentum is defined to be a momentum at which the energy equals to 
the chemical potential at zero temperature: /i = e{pp), and the Fermi velocity is defined as 



dtip) 

vf = 



(10) 



P=Pf 



<^ Pf, we may integrate out the fast modes to get an 



If we are interested in low energies, 
effective operator, 

Meft = -dr -VF-l, (11) 

which has complex eigenvalues. However, when we include the —vp sector, we have 
Meft{vf)Meft{—'Vf) < (i.e., has real negative eigenvalues), assuming e{p) = e{—p). 
We again see that the sign problem in the electron system is alleviated in the low-energy 
effective theory. 

III. QCD 

Let us recall why the measure of dense QCD is complex in Euclidean space. We use the 
following analytic continuation of the Dirac Lagrangian to Euclidean space: 

xo -ix%, Xi Xe ; 7o If, li ^ ilE ■ (12) 
The Euclidean gamma matrices satisfy 

rJ = l'E , {7^^,7^1 = 25'^^ (13) 
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The Dirac-conjugated field, 'ip = V'W^; is mapped into a field, still denoted as t/S, which is 
independent of ip and transforms as ip'^ under SO {4). Then, the grand canonical partition 
function for QCD is 

Z{^j)= j dA^Aei{M)e~^^^^\ (14) 

where S{A^) is the positive semi-definite gauge action, and the Dirac operator 

M = 7^D^ + /x7^, (15) 

where De = dE + i^E is the analytic continuation of the covariant derivative. The Hermitian 
conjugate of the Dirac operator is 

Mt = -7^D^ + ;x7^ . (16) 

The first term in ()15|) is anti-Hermitian, while the second is Hermitian, hence the generally 
complex eigenvalues. When /z = 0, the eigenvalues are purely imaginary, but come in 
conjugate pairs (A, A*) 11], so the resulting determinant is real and positive semi-definite: 

det M = JJ A*A > . (17) 

In what follows we investigate the positivity properties of an effective theory describing 
only modes near the Fermi surface. A system of degenerate quarks with a net baryon number 
asymmetry is described by the QCD Lagrangian density with a chemical potential /i, 

/:qcd = i^ilpi^- \f;,F'^^' + /x^7oV', (18) 

where the covariant derivative = + and we neglect the quark mass for simplicity. 

At high density (/x 3> Aqcd), due to asymptotic freedom the energy spectrum of quarks 
near the Fermi surface is approximately given by a free Dirac eigenvalue equation, 

(a ■p-/i)?/'± = E±?/)±, (19) 

where a = 707 and 'ip± denote the energy eigenfunctions with eigenvalues E± = — yU ± 
respectively. At low energy E < fi, the states 1^+ near the Fermi surface, \p\ ~ /i, are easily 
excited but ip-, which correspond to the states in the Dirac sea, are completely decoupled 
due to the presence of the energy gap /i provided by the Fermi sea. Therefore the appropriate 
degrees of freedom at low energy consist of gluons and only. 

Now, we wish to construct an effective theory describing the dynamics of by integrating 
out modes whose energy is greater than /i. Consider a quark near the Fermi surface, whose 
momentum is close to ^vp- Without loss of generality, we may decompose the momentum 
of a quark into a Fermi momentum and a residual momentum as 

Pi, = fJ-v^ + (20) 
where v'^ = (0, vp)- Since the quark energy is given as 

E = -f,+ y/(/||+/i)2 + /i, (21) 
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the residual momentum should satisfy (/y + yu)^ + l± < 4:fi^ with /y = vpl ■ vp and l± = I — 

To describe the small excitations of the quark with Fermi momentum, fiVp, we decompose 
the quark fields as 

^(x) = e'^""^-^ [iP+{vF, x) + ^.{VF, x)] , (22) 

where 

iJi{vp,x) = P^{vp)e-'^'^--'ij{x) with P^(vp)^^^^^. (23) 
The quark Lagrangian in Eq. (fTK|) then becomes 

ijj {ip+ /i7°) ip = i>+{vF, x)i'yil^D^^/j+{vF, x) + i'-ivF, a;)7° (2/i + iDy) ip-^{vF, x) 

+ [tp-ivF, x)ip^ip+{vF, x) + h.c] (24) 

where 7f = (7°, vfVf ■ 7), 7± = 7^ - 7f , ^|| = Vf'D^ with V = (1, tT^), V"'' = (1, -vf), and 

At low energy, we integrate out all the "fast" modes ip- and derive the low energy effective 
Lagrangian by matching all the one-light-particle irreducible amplitudes containing gluons 
and in loop expansion. The effects of fast modes will appear in the quantum corrections 
to the couplings of low energy interactions. At tree-level, the matching is equivalent to 
eliminating in terms of equations of motion: 

Therefore, the tree-level Lagrangian for ip^ becomes 

= V^+^7p^^+ - ^^+7°(^x)'V^+ + • • • , (26) 

where the ellipsis denotes terms with higher derivatives. 

Consider the first term in our effective Lagrangian, which when continued to Euclidean 
space yields the operator 

Meft=7f-^(^). (27) 

Meft is anti-Hermitian and it anti-commutes with 75, so it leads to a positive semi-definite 
determinant. However, note that the Dirac operator is not well defined in the space of 
ip+{vF,x) (for fixed vf), since it maps ip+i'VpyX) into ip+i~'^F,x): 

ilpW P+i) = P-iIp\\ ifj. (28) 

Since P_(t7p) = P+{—vf), ip ip+i^vpyx) are 4'+{~'Vf,x) modes, or fluctuations of a quark 
with momentum —fivp- 

We can demonstrate the necessity of including both ip+{vF,x) and ijj+{—VF,x) modes in 
our effective theory by considering charge conservation in a world with only +Vf quarks. 
The divergence of the quark current at one loop is at the leading order in the 1/fi expansion 

{d,r^'{vF, x)) =gsj ^e-^- Vn J(p)A;^(-p) , (29) 
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where A\\ = {Aq, vpVp ■ A) and 11^^ is the vacuum polarization tensor in the effective theory 
given as 0] 

Kb (P) = --^^^ab — . . . . ^ +91 ] , 30) 

where g'^ = g'^'^ — {V'^V'^ + V^V'^)/2 is the metric tensor perpendicular to and V'^. The 
polarization tensor has to be transverse to maintain gauge invariance. We find that if we 
have both fields ipj^{vF, x) and ?/'+(— -uf, x) the current is conserved and the gauge symmetry 
is not anomalous: 

{d,Jli{vF, x) + d.JHi-VF, x)) = 0. (31) 

Therefore, we need to introduce quark fields with opposite momenta. The Dirac operator is 
well defined on this larger space. 

This anomaly can be understood in terms of spectral flow, since the Fermi surface is 
(in a certain sense) not gauge- invariant. Under a gauge transformation, U{x) = e*'''^, the 
Hamiltonian changes and the energy spectrum of free modes of residual momentum / shifts 
to E = I ■ vp + q ■ vf- Quarks near the Fermi surface with vf ■ q> flow out of the Fermi 
sea, creating charge. This charge creation is compensated by quarks with opposite vp] their 
energy decreases and they flow into the Fermi sea. However, unless modes with opposite 
velocities (i.e. both sides of the Fermi sphere) are included, charge is not conserved. 

Thus far we have considered the quark velocity as a parameter labelling different sectors of 
the quark field. This is similar to the approach of heavy quark effective theory (HQET) , 
in which the velocity of the heavy charm or bottom quark is almost conserved due to 
the hierarchy of scales between the heavy quark mass and the QCD scale. However, this 
approach contains an ambiguity often referred to as "reparameterization invariance" , related 
to the non-uniqueness of the decomposition ()20|) of quark momenta into a large and residual 
component. In the dense QCD case, two ip{vF,x) modes whose values of vp are not very 
different may actually represent the same degrees of freedom of the original quark field. In 
what follows we give a different formulation which describes all velocity modes of the quark 
field, and is suitable for defining the quasiparticle determinant. 

First, a more precise definition of the breakup of the quark field into Fermi surface modes. 
Using the momentum operator in a position eigenstate basis: p = —id, we construct the 
Fermi velocity operator: 

^ —i d , , 

which is Hermitian, and a unit vector. 

Using the velocity operator, we define the projection operators P± as before and break 
up the quark field as, ipi^x) = ip+^x) + ip-{x), with ip± = P±ip. By leaving v as an operator 
we can work in coordinate space without introducing the HQET-inspired velocity Fourier 
transform which introduces vp a.s a parameter. If we expand the quark field in the eigenstates 
of the velocity operators, we recover the previous formalism with all Fermi velocities summed 
up. 

The leading low-energy part of the quark action is given by 

C+=^P^{v){z^-4 + fi^o)P+{v)^ . (33) 
As before, we define the fields ip-^ to absorb the large Fermi momentum: 

^+{x) = e-''''^-'^P+{v)tlj{x). (34) 
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Let us denote the eigenvalue v obtained by acting on the field (which has momentum 
of order jj) as vi (or v "large"), whereas eigenvalues obtained by acting on the effective 
field theory modes ■?/'+ are denoted Vr (or v "residual"). If the original quark mode had 
momentum p with \p\ > jj, (i.e. was a particle), then vi and Vr are parallel, whereas if \p\ < fj, 
(as for a hole) then Vr and vi are anti-parallel. In the first case, we have P+{vi) = P^{vr) 
whereas in the second case P+{vi) = P^{vr)- Thus, the residual modes can satisfy either 
of P±{vr)'4^+ = depending on whether the original ip mode from which it was derived 
was a particle or a hole. In fact, ip^ modes can also satisfy either of P±{vi)tp+ = tp^ since 
they can originate from ip modes with momentum ~ +fiv as well as —fiv (both are present 
in the original measure: Dip Dip). So, the functional measure for ?/^+ modes contains all 
possible spinor functions - the only restriction is on the momenta: |/o|, |^| < A, where A is 
the cutoff. 

In light of the ambiguity between vi and w^, the equation ip = e~^^^^'^ip+ must be modified 

to 

ip = exp {+ifix ■ V a ■ v) tp+ = exp {+ifix ■ Vr a ■ Vr) ip+ , (35) 

where the factor of a ■ Vr corrects the sign in the momentum shift if Vr and vi are anti- 
parallel. In general, any expression with two powers of v is unaffected by this ambiguity. 
For notational simplicity we define a local operator 

X ^ ^x-va-v = l^^g-g-- (36) 

Taking this into account, we obtain the following action: 

/:+ = V^+e-*^(z^-4 + /i7o)e+^^^+ . (37) 

We treat the ^ term separately from + /i7o since the former does not commute with X, 
while the latter does. Continuing to Euclidean space, and using the identity P_7^P+ = 7II-P+, 
we obtain 

= {d^ + tA^^) ^+ , (38) 

where 

A>1 = e-'^ A^ e+'^ , (39) 

and all 7 matrices are Euclidean. The term containing A cannot be fully simplified because 
[v, A] 7^ 0. Physically, this is because the gauge field carries momentum and can defiect the 
quark velocity. The redefined modes are functions only of the residual momenta 1, and 
the exponential factors in the A term refiect the fact that the gluon originally couples to the 
quark field not the residual mode ip+- 

The kinetic term in ()38|) can be simplified to 

= (40) 

since v ■ dv ■ 'y = (9-7 . The action (|HHj) is the most general dimension 4 term with the 
rotational, gauge invar iance j3| and projection properties appropriate to quark quasipar- 
ticles. Therefore, it is a general consequence of any Fermi liquid description of quark-like 
excitations. 

The operator in (|38|) is anti-Hermitian and leads to a positive, semi-definite determinant 
since it anti-commutes with 75. The corrections given in (|26p are all Hermitian, so higher 
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orders in the expansion may re-introduce complexity. The structure of the leading 
term plus corrections is anti-Hermitian plus Hermitian, just as in the original QCD Dirac 
Lagrangian with chemical potential. 

By integrating out the fast modes, the Euclidean QCD partition function can be rewritten 

as 

Z{fx) = j dA+ detMefT(A+)e-^°ff(^+\ (41) 

The leading terms in the effective action for gluons (these terms are generated when we 
match our effective theory, with energy cutoff A, to QCD) also contribute only real, positive 
terms to the partition function: 

S,^{A) = I d'xE (If^.F^^ + ^ E ^1.^1 J > 0' (42) 

\ Vp J 

where = A — A\\ and the Debye screening mass is M = N f / {2'n'^)gs^ . Note that 
Landau damping is due to softer quark modes which have not been integrated out, and 
therefore do not contribute to matching. 

Although the HDET only describes low-energy modes, it still contains Cooper pairing 
interactions. This is because Cooper pairing, in which the quasiparticles have nearly equal 
and opposite momenta, is induced by gluonic interactions with small energy and momentum 
transfer. That is, although a gluon exchange (or other interaction) which causes a large 
angular deflection of a quasiparticle 

\P) - \P') 

must involve a large momentum transfer, and hence is not part of the effective theory, a 
Cooper pairing interaction 

|p, -p) \p\ -p) 

only involves a small energy and momentum transfer, even if the angle between p and p' is 
large. Hence, it is described by the leading order interaction between soft gluons and quarks 
in the effective theory (jHHj) . 

Matching of hard gluon effects also leads to four-quark operators in the effective theory. 
The addition of these four-quark operators still leads to a positive action for attractive chan- 
nels, since they arise from quasiparticle-gluon interactions which are originally positive. A 
simple way to study the positivity of four-quark operators is to replace them by a vector 
field with trivial quadratic term V'^ which couples to quarks like the original gluon: V^ip'^^il). 
Completing the square, we see that the resulting path integral is positive. Note that this 
argument does not apply to interactions involving six or more quarks, or interactions in- 
volving virtual anti-quarks. However, these are always suppressed by additional powers of 
/i. 

Finally, some comments on more general models than QCD - in particular electronic 
systems in condensed matter - which are also affected by a sign problem. We note that as 
long as the important dynamics of the system involve energies which are small compared to 
the Fermi energy (or are exhibited at a temperature small compared to the Fermi temper- 
ature), the Euclidean description of the system is likely to have a positive effective action 
after modes far from the Fermi surface are integrated out. Models involving, e.g., attractive 
four-fermion interactions or long range gauge fields fall into this category. As long as the 
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important interactions are soft interactions between quasiparticle degrees of freedom not too 
far from the Fermi surface, the model can be simulated without a sign problem, even if it is 
strongly coupled. 

As a specific example, consider the Hamiltonian 

H = tabc\cb + ^ Vabcd c\ci,clcd , (43) 

ab abed 

where abed denote lattice site and spin indices. This Hamiltonian includes the Hubbard 
model as a special case. A constant diagonal part of the kinetic term acts as a chemical 
potential yU = taa-, but t may also contain hopping terms such as a nearest-neighbor off- 
diagonal term in the tight binding approximation. The usual kinetic term, taken by itself, 
produces a band structure which in momentum space has the form 

Y,E{k)clck , (44) 

k 

for some function E{k). A non-zero chemical potential leads to filling of levels up to some 
momentum kp. There is an effective theory description of modes near this surface if the 
Fourier transform of the interaction term Vabcd only contains soft interactions (i.e. support 
for typical momentum transfers less than kp.) 

To obtain the partition function Z we must extend the fields from the d-dimensional 
lattice to (d+l)-dimensional Euclidean space, with finite extent (3 in the time direction. 
The action density is given by 

(«) 

where Ti is the Hamiltonian density. Using the Hubbard-Stratonovich transformation, we 
can rewrite the interaction term in terms of a functional integral over a background field U. 

Z = ^ JYI dUab e^^^^f ^^"''c.lfUef det [ T - f/ ] , (46) 

ab 

where T = 4- — t contains the time derivative as well as the kinetic term t from the 



Hamiltonian ()43|) . We assume that the interaction v is real and negative definite so that the 
exponential is real and the integral converges. 

The determinant can be expanded about our effective theory. The leading term is positive 
and has no sign problem. The correction terms are suppressed as long as the background 
field U has support on momentum scales small compared to kp. (See discussion in next 
section, especially equation (fS^ with U playing the role of the gauge field). Typical fields 
U in the integral are determined by the exponential term - if the Fourier transform of v has 
little support at large momenta, then the exponent U{k) v~^{k) U{k) will be large and 
negative, and the corresponding modes of U suppressed. 

Unfortunately, the Hubbard model itself contains the interaction term 



Ui^rii^ (47) 



where nj = c|cj is the number operator. This interaction allows all momentum transfers 
up to Tc/a, where a is the lattice spacing. Because it couples low-energy modes to fast 
modes far from the Fermi surface, the resulting partition function isn't well-approximated 
by the positive EFT part. However, a model with modified (softer) interactions would be 
approximately positive and might in fact be more physically realistic. 
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IV. LATTICE SIMULATION 



The goal of this section is to give a method for simulating QCD at finite density. We 
will consider a chemical potential fi much larger than Aqcd throughout, and divide the 
functional integral over quark excitations into two parts: (I) modes within a shell of width 
A of the Fermi surface, and (II) modes which are further than A from the Fermi surface. 
We will assume the hierarchy 

;u » A » Aqcd • (48) 

The quark determinant in region (I) is well approximated by the determinant of the 
leading operator in high density effective theory (HDET) as long as the first inequality in 
PH|) is satisfied. As discussed in the previous section, it is positive and real. 

Here we will show that the contributions to the effective action for the gauge field from 
quark modes in region (II) are small and vanish as the A grows large compared to Aqcd- 

First consider the theory in Minkowski space. The Dirac operator is 

M = ip+ fi-fo (49) 
and the Dirac equation can be written as 

idotp = Hil) (50) 

with 

H = ia ■ d — (51) 

a Hermitian operator. The break up into regions (I) and (II) proceeds naturally in terms 
of energy eigenvalues of H (or Iq in the HDET notation). The low- lying modes in region (I) 
are particle states with spatial momenta satisfying \p\ ^ /x. 

The analytic continuation of region (I) to Euclidean space leads to the HDET determinant 
considered previously. 

Modes in region (II) all have large energy eigenvalues, at least as large as A. In considering 
their effect on physics at the scale Aqcd, we can integrate them out in favor of local operators 
suppressed by powers of Aqcd/A. 

To make this concrete, consider the effective action for gauge fields with field strengths 
F^i, of order Aqcd- The quark contribution to this effective action is simply the logarithm of 
the determinant we wish to compute. It can be expanded diagrammatically in graphs with 
external gauge field lines connected to a single quark loop. Restricting to region (II), we 
require that the quark modes in the loop have large H eigenvalues. Evaluating such graphs 
leads only to operators which are local in the external fields A^(x). 

The resulting renormalizable (dimension 4) operator is the finite density equivalent of F^^^, 
except that due to the breaking of Lorentz invariance it contains separate time- and space-like 
components which represent the contribution of high-energy modes to the renormalization 
of the coupling constant, and Debye screening. These effects do not introduce a complex 
component when continued to Euclidean space. 

Higher dimension operators, which involve additional powers of F^,^ or covariant deriva- 
tives Dfj^ are suppressed by the scale A. These are presumably the source of complex terms 
introduced to the effective action. However, due to the 1/A suppression they are dominated 
by the contribution from the low- lying modes in region (I), which is necessarily non-local, 
but real. 
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The logarithm of the Euchdean quark determinant will have the form: 

IndetM ~ + (non — local, real) + (local, complex) , (52) 

where the first term is the (real, constant) free energy of non-interacting quarks, the second 
term is from the positive determinant in region (1) and the last term is the suppressed, 
complex contribution from region (II). Only the last two terms depend on the gauge field 

On the lattice, one can use the dominant dependence of detM on the first and second 
terms to do importance sampling. In order to keep the complex higher dimension operators 
(last term in ()52|) ) small, it is important that the gauge field strengths are kept smaller 
than A^. One can impose this condition by using two different lattice spacings, Qg for the 
gluons and Odet for the quarks, with ag > Odet- The determinant is calculated on the finer 
Odet lattice, and is a function of plaquettes which are obtained by interpolation from the 
plaquettes on the coarser ag lattice. Interpolation can be defined in a natural way, since 
each lattice link variable U^fj, is an element of the gauge group, and one can connect any two 
points gi, g2 on the group manifold in a linear fashion: g{t) = gi+ t{g2 — gi) , <t < 1. 

More explicitly, let x and x + ttgfi be two neighboring points on the coarse lattice, and 

Zn = X + ncdet/^ , n = 0,1, ■ ■ ■ N (53) 

be the corresponding points on the fine lattice: zq = x, zn = x + ttgfi, where N = Og/adet- 
Then a link 

U,,^ = exp[tagG'^n 

is interpolated to a set of links as 

=exp[mdetm"^] , (54) 

where the t"* are SU(N) generators, and the bar denotes the finer lattice. Equation (j54p 
allows us to compute f7 links which are sub-links of original ag link variables, and lie on the 
outer perimeter of a plaquette. The remaining U links, which are within an ag plaquette, can 
be obtained through a similar interpolation starting from opposite sides of the perimeter, 
yielding an entire set of plaquettes on the Odet lattice. The field strengths resulting from 
this interpolation are always of order and can be kept small compared to the cutoff 
A^. To properly include the quasiparticle modes, the spacing of the fermion lattice must 
be Odet ~ l/yW, while ag ~ 1/Aqcd is probably sufficient to capture the effects of non- 
perturbative gauge configurations. 

The fermion determinant is to be computed as a function of the finer plaquettes {Ifxti}- 
The result is (approximately) real and positive and can be used for importance sampling. 
Further, there is a physical understanding of the complex part of the determinant: it origi- 
nates in the modes far from the Fermi surface which have been integrated out. 

V. INEQUALITIES AND ANOMALY MATCHING 

Positivity of the measure allows for rigorous QCD inequalities at asymptotic density. For 
example, inequalities among masses of bound states can be obtained using bounds on bare 
quasiparticle propagators. One subtlety that arises is that a quark mass term does not lead to 
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a quasiparticle gap (the mass term just shifts the Fermi surface). Hence, for technical reasons 
the proof of non-breaking of vector symmetries jl^] must be modified. (Naive apphcation 
of the Vafa-Witten theorem would preclude the breaking of baryon number that is observed 
in the color-fiavor-locked (CFL) phase [l^)- A quasiparticle gap can be inserted by hand 
to regulate the bare propagator, but it will explicitly violate baryon number. However, 
following the logic of the Vafa-Witten proof, any symmetries which are preserved by the 
regulator gap cannot be broken spontaneously. One can, for example, still conclude that 
isospin symmetry is never spontaneously broken (although see below for a related subtlety). 
In the case of three flavors, one can introduce a regulator d with the color and flavor structure 
of the CFL gap to show rigorously that none of the symmetries of the CFL phase are broken 
at asymptotic density. On the other hand, by applying anomaly matching conditions jl6| . 
we can prove that the SU{?>)a symmetries are broken. We therefore conclude that the CFL 
phase is the true ground state for three light flavors at asymptotic density, a result that was 
first established by explicit calculation [oT 17. isl . 

To examine the long-distance behavior of the vector current, we note that the correlator 
of the vector current for a given gauge field A can be written as 

{j;{vF,x)Jl{vF,y))^ = -Tr7^T^^^(x,y;c?)7.r''S^(z/,x;rf), 

where the SU{Nf) fiavor current J^{vf, x) = ip+{vF, x)'y^T""ilj+{vF, x). The propagator with 
S'[/(3)\/-invariant IR regulator d is given as 

S^x, y- d) = {x\ - \y) = dr (x| e^^^-*^) \y) 



where with D = d + iA 



M = Jo 



D-V d _ 
d^ D-Vj- 



Since the eigenvalues of M are bounded from below by rf, we have 



M 



-dR 



■^{x\x)^{y\y), 



(55) 



where R= \x ~ y\. The current correlators fall off rapidly as -R — > oo; 



dA+ detMeff(A) e-^=ff {j^{vF,x)J^{vF,y)y 



< 



'X{vF,x)J^{vF,y)y 



< 



-2dR 



d^ 



\{x\x)\ \{y\y)\ 



(56) 



where we used the Schwartz inequality in the first inequality, since the measure of the ef- 
fective theory is now positive, and equation in the second inequality. The IR regulated 
vector currents do not create massless modes out of the vacuum or Fermi sea, which implies 
that there is no Nambu-Goldstone mode in the SU{3)v channel. Therefore, for three mass- 
less fiavors SU{3)v has to be unbroken as in CFL. The rigorous result provides a non-trivial 
check on explicit calculations, and applies to any system in which the quasiparticle dynamics 
have positive measure. 
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It is important to note the order of limits necessary to obtain the above results. Because 
there are higher-order corrections to the HDET, suppressed by powers of A//i, that spoil its 
positivity, there may be contributions on the RHS of of the form 

O (^) f{R) , (57) 

where f{R) falls off more slowly than the exponential in (j56|) . To obtain the desired result, 
we must first take the limit fi oo at fixed A before taking R oo. Therefore, our results 
only apply in the limit of asymptotic density. 

Although our result precludes breaking of vector symmetries at asymptotic density in the 
case of three exactly massless quarks j2l| , it does not necessarily apply to the case when the 
quark masses are allowed to be slightly non-zero. In that case the results depend on precisely 
how the limits of zero quark masses and asymptotic density are taken, as we discuss below. 

In fl3\ the authors investigate the effect of quark masses on the CFL phase. These 
calculations are done in the asymptotic limit, and are reliable for sufficiently small quark 
masses. When = m^^ = m << (unbroken SU{2) isospin, but explicitly broken 
SU{3)), one finds a kaon condensate. The critical value of at which the condensate 
forms is m* ~ m ^/3Ao/^ where Aq is the CFL gap (see, in particular, equation (8) of the 
first paper). As kaons transform as a doublet under isospin, the vector SU{2) symmetry is 
broken in seeming contradiction with our result. 

However, a subtle order of limits is at work here. For simplicity, let us set m = 0. Note 
that the CFL regulator d, which was inserted by hand, explicitly breaks SU{3)a through 
color-fiavor locking, leading to small positive mass squared for the pions and kaons, given as 



,2 . . „, j2 



a J' In . (58) 

The meson mass is not suppressed by since, unlike the Dirac mass term, the regulator, 
being a Majorana mass, does not involve antiquarks 0]. 

Therefore, even when the light quarks are massless, there is a critical value of rris necessary 
to drive negative the mass-squared of kaons and cause condensation: 



"^7,rf/iln(0]'^' > (rf/i)V2, (59) 



where Qs is the strong coupling constant. Note the product of Qs with the logarithm grows 
as n gets large. To obtain our inequality we must keep the regulator d non-zero until the 
end of the calculation in order to see the exponential fall off. To find the phase with kaon 
condensation identified in we must keep rris larger than m*. (Note /i oo, so to have 
any chance of finding this phase we must take d —>■ keeping dR large and dfi small.) 
Since the UV cutoff of the HDET must be larger than m^, we have 

l>(f)^l^ (ao) 

which implies 

^f{R) > ^f{R). (61) 
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Note the right hand side of this inequahty does not necessarily fall off at large R, and also 
does not go to zero for /i — oo at fixed A and d. This is a problem since to apply our 
inequality the exponential falloff from ()56p must dominate the correction term (|37j). which 
is just the left hand side of (|HT|). Combining these equations, we see that the exponential 
falloff of the correlator is bounded below, 

-2dR J 

— >jm- (62) 

in the scaling region with a kaon condensate, > m*. 

Alternatively, if we had taken m<j to be finite for fixed regulator d (so that, as /i — ^ oo, 
eventually < m*), the inequality in could be applied to exclude a Nambu-Goldstone 
boson, but we would find ourselves in the phase without a kaon condensate. 



VI. CONCLUSION 

The low-energy physics of dense fermionic matter, ranging from quark matter to electronic 
systems, is controlled by modes near the Fermi surface. An effective Lagrangian describing 
the low-energy modes can be given in a systematic expansion in powers of the energy scale 
over the chemical potential. The leading term in this expansion has a simple form, and we 
have shown that it leads to a real, positive Euclidean path integral measure. 

This observation opens the door to importance sampling in Monte Carlo simulations of 
dense matter systems. The key requirement is that the interactions do not strongly couple 
the low-energy modes to modes far from the Fermi surface. QCD at high density satisfies this 
requirement, as do all asymptotically free models. Electronic system in which the important 
interactions involve momentum transfer less than the Fermi energy are in this category, 
although some idealized models such as the Hubbard model are not. We have given some 
proposals for how the positive effective theory might be simulated numerically. Ultimately, 
we hope that actual practitioners will develop even more practical methods. 

Finally, positivity has analytical applications as well, since it allows the use of rigor- 
ous inequalities. In QCD we obtain restrictions on symmetry breaking patterns at high 
density. Similar restrictions can probably be obtained for electronic systems with suitable 
interactions. 
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